Cassini UVIS Observations of the lo Plasma Torus. 
IV. Modeling Temporal and Azimuthal Variability 



A. J. Steffi 

Southwest Research Institute, Department of Space Studies, 1050 Walnut St., Suite 300, 

Boulder, CO, 80302 

stef f KSboulder . swri . edu 

P. A. Delamere, F. Bagenal 
University of Colorado, Laboratory for Atmospheric and Space Physics, Campus Box 392, 

Boulder, CO 80309-0392, USA. 

Received ; accepted 



- 2 - 
ABSTRACT 

In this fourth paper in a series, we present a model of the remarkable tempo- 
ral and azimuthal variability of the lo plasma torus observed during the Cassini 
encounter with Jupiter. Over a period of three months, the Cassini Ultravio- 
let Imaging Spectrograph (UVIS) observed a dramatic variation in the average 
torus composition. Superimposed on this long-term variation, is a 10.07-hour 
periodicity caused by an azimuthal variation in plasma composition subcorotat- 
ing relative to System III longitude. Quite surprisingly, the amplitude of the 
azimuthal variation appears to be modulated at the beat frequency between the 
System III period and the observed 10.07-hour period. Previously, we have suc- 
cessfully modeled the months-long compositional change by supposing a factor of 
three increase in the amount of material supplied to lo's extended neutral clouds. 
Here, we extend our torus chemistry model to include an azimuthal dimension. 
We postulate the existence of two azimuthal variations in the number of super- 
thermal electrons in the torus: a primary variation that subcorotates with a 
period of 10.07 hours and a secondary variation that remains fixed in System III 
longitude. Using these two hot electron variations, our model can reproduce the 
observed temporal and azimuthal variations observed by Cassini UVIS. 

Subject headings: Jupiter, Magnetosphere; lo; Ultraviolet Observations; Spectroscopy 
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Introduction 



During the Cassini spacecraft's flyby of Jupiter (October 2000 through March 2001^ 
the Ultraviolet Imaging Spectrograph (UVIS) made extensive observations of the lo 



plasma torus. 



(lEsposito et al. 



he s ensitivity, bandpass, resolution, and imaging capabilities of UVIS 



2004 ) coupled with the temporal coverage of the observations make this a 
particularly rich dataset. Analysis and modeling of these observations has led to remarkable 
new insights into the behavior of the lo torus. In this paper we present the results of our 
efforts to model the ter nporal and azimut hal variability of the lo plasma torus observed by 



UVIS and discussed by 



Steffletal 



(120061 ). To put this effort into proper context, here we 



recapitulate the prior analysis and modeling of the Cassini UVIS observations. 

Information about the UVIS lo torus dataset, including examples of the observing 
geometry, images of the raw and processed data, and des criptions of the da ta reduction 



and calibration techniques used was presented in paper I (ISteffl et al. 



2004al ). The analysis 



of a small subset of UVIS lo torus observations obtained in January 2001, sho rtly after 



the Ca ssini spacecraft's closest approach to Jupiter, w as presented in Paper II (ISteff 



2004bl ). Using the CHIANTI atomic physics database (jPere et al. 



1997 



Young et al 



et al 



20031) 



they derived the composition and electron temperature of torus plasma as a function 
of radial distance from Jupiter. The torus composition presented in Pap er II differs 



significantly from that derived from observations made in the Voyager ei a. ( lBagenallll994l ). 
with less oxygen and a lower electron temperature. The CHIANTI-based spectral model of 
the lo torus was used to analyze UVIS data obtained during a 45-day period on the Cassini 
spacecraft's approach to Jupiter. Significant temporal and azimuthal variati ons in torus 



composition were found, which were presented in paper III (jSteffl et al. 



20061). 



To understand the processes governing the lo plasrna toru s and their changes between 



the Voyager and Cassini epochs iDelamere and Bagena] 



(I2OO3I ) (hereafter referred to as 
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DB03) developed a physical chemistry model of the 



mode 



ing w ork by 



(119981), and 



Barbosa et al. 



Lichtenberg et al. 



19831), 



ShemanskyI (119881 ) . 



o torus which b uilds o n previous 



Barbosal fll994h . 



Schreier et al 



(l200ll ). The DB03 model, from which the azimuthal models 



discussed below are directly descended, is a "0-d" or one-box model that calculates the 
flow of mass and energy through one cubic centimeter of torus plasma placed at a radial 
distance of 6 Rj. The effects of electron impact ionization, recombination (both radiative 
and dielectronic), charge exchange reactions, parameterized radial transport. Coulomb 
collisions between species, and radiative energy losses on the torus plasma are all included. 

New mass is supplied to the torus in the form of neutral oxygen and sulfur atoms 
which eventually become ions through either electron impact ionization or charge exchange 
reactions. Conversely, mass is lost from the model when an ion becomes neutralized 
through charge exchange or recombination or through outward radial transport. The 
rotation speed of plasma in the lo torus (~75 km/s at 6 Rj) significantly exceeds Jupiter's 
escape velocity. When torus ions becomes neutralized, the resulting atom is no longer 
constrained by Jupiter's magnetic field and is quickly lost fr om the lo torus, even tually 



forming an extended nebula hundreds of Jovian radii in size (IMendillo et al. 



2004h . Torus 



plasm a is also convected radially outwards by the interchange motions of magnetic flux 



tubes ( iRichardson and Siscoe 



1981 



Siscoe and Summers 



198ll ). However, the details of 



convective radial transport are beyond the scope of the DB03 model; instead, the lifetime 
of torus plasma against radial transport is specified by an input parameter, tq. 

Energy is supplied to the torus via the "pickup energy" imparted to new ions as they 
are accelerated fro m Keplerian orbit a l velocities to near co-rotational with the magnetic 



field. As noted by 



ShemanskyI p98 



Smith et al. 



(119881 ). and others, pickup energy alone 
cannot supply the torus with enough energy to maintain roughly 10^^ W of UV radiation. 



a thermal electron temperature of ~5 eV, and an average ionization state of ~1.5, as are 
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observed: an additional source of energy is required. It is now widely accepted that this 
additional energy is provided by a small population of super-thermal electrons. This is 
supported by detections of a h igh-energy component of the torus electron d istribution 



function by both the 



(iFrank and Paterson 



Voya ger ( IScudder et al. 



1981 



Sittler and Strobel 



19871 ) and Galileo 



20001 ) plasma instruments. In addition, analysis of Ulysses URAP 



observations showe d the electron distribution function in the lo plasma torus resembles a 



kappa distribution flMeyer-Vernet et al 



19951) 



In the torus model of DB03, the electron distribution function is simplified as the 
sum of two Maxwellian populations: a thermal population near the canonical 5 eV and 
a small (roughly 0.2% of the total electron density), hot pop ulation with temperature 



between 50-100 eV, consistent with 



Sittler and Strobell (119871 ). The temperature of the hot 



population is held constant at a value specified by the input parameter, T^^hot- Since the 
hot electron population rapidly couples to the thermal electron population via Coulomb 
collisions (with a characteristic cooling time on the order of 30 minutes), maintaining a 
constant hot electron temperature requires the hot population be continuously resupplied 
with energy. For the plasma conditions observed by UVIS during the Cassini epoch, DB03 
found the hot electron source was responsible for up to 60% (10^^ W) of the total energy 
input to the torus. 

The energy for these hot electrons must ultimately be derived from Jupiter' s rotation , 



but th e particular details of the heating mechanism remain poorly understood. 



Barbosa 



(119851 ) suggested that the hot electrons may be heated by lower hybrid waves generated 



during the thermalization of the pickup io n ring beam distribution. However, if lo's extende d 



neutral clouds are highly peake d near lo (jSmyth and Marconi 



2003 



2004 



Smyth and Marconi 



Burger and Johnson 



20051 ). this mechanism will produce a 13- hour periodicity (the 



synodic period between lo and System III) which should be detectable in the Cassini 
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UVIS observations. No such 13-hour periodicity was seen, suggesting that rin 



Thorne et al 



g beam 



(I1997I ) reported 



thermahzation is not the dominant production mechanism, 
rapid (100 km/s) inward motion of hot tenuous plasma during a flux tube interchange 
event. Such rapid motion could supply the lo torus with electrons that have been heated 
in the outer torus/middle magnetosphere. Alternatively, we propose that the hot electrons 
are produced locally throughout the torus during small-scale flux tube interchange events. 
These interch ange events likely generat e the field-aligned electron beams detected by the 



Galileo PLS (IFrank and Paterson 



2OOOD in the shear region between inward and outward 



moving flux tubes. 

An estimate of the total power available to the torus from these beams of field-aligned 
electrons can be made using the following equation: 



where $ is the flux of electrons in the field-aligned beams, T is the average electron energy, 
/ is the fraction of the torus in which the be ams occur, and r^ and a re radial distance of 



the inner and outer edges of the torus. From 



Frank and PatersonI (120001 ) . $ ^ 10*^ cm"^ s~^ 



T 500 eV (ranging between 100 eV and a few keV), and / ~ 0.1. Assuming a torus that 
extends from 6.0 to 7.5 Rj, yields a total power input of a few 10^^ W, consistent with what 
is required by DB03. 

While a discussion of the origin of these beams is beyond the scope of this paper, we 
note that a plausible heating mechanism is the propaga tion of inertial Alfven waves out of 



the plasma torus, as was first discussed by 



Craryl (119971 ) for the lo fl ux tube. We note tha t 



the en ergy spectrum (100 eV to few keV) of the beams reported by iFrank and Paterson 



(120001 ) is consistent with electron acceleration to roughly the Alfven velocity just outside of 



the torus. For additional discussion of acceleration by inertial Alfven waves see 



Su et al. 
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(120061) and lSwiftI (120061 ). 



1.1. Temporal variability during the Cassini era 



Analysis of the Cassini UVIS observations of the lo torus showed signific ant changes 



in th e composition of torus plasma between October 2000 and January 2001 (IStefB et al. 



20061 ) ■ During this period, the average mixing ratio (ion density divided by electron density) 
of S II in the torus declined by a factor of 2 during this period, with a corresponding factor 
of 2 increase in the mixing ratio of S IV. Observations by the Galileo Dust Detector System 
(DDS) showed a dramatic increase, by over three orders of magr iitude, in the emiss ion 



20031). The 



rate of logenic dust immediately prior to the UVIS observations (iKriiger et al. 
enhanced dust emissions began in July 2000, peaked in September 2000, and returned to 
pre-event levels by December 2000. Presumably, this increase (the largest, by an order 
of magnitude, observed by the Galileo DDS over a seven-year period) was in response to 
some major volcanic event on lo, possibly the eruption of Tvashtar Catena that deposited 
a Pele-like ring of red ma terial 1200 km in d iameter on lo's surface between February 



2000 and December 2000 flGeissler et al 



20041 ). A 400 km-high plume o ver this region wa s 



observed by the Gassini Imaging Science Subsystem in December 2000 (jPorco et al. 



20031). 



These observations led 



Delamere et al 



(120041 ) to propose that an increase in the amount 



of material supplied to the lo neutral clouds might be responsible for the compositional 
changes seen by Cassini UVIS. Their model included a Gaussian increase in the neutral 
source rate with an amplitude of 3.5 and a width of 22.5 days. As the supply of neutrals 
increased, the densities in the model torus rapidly increased to unrealistic levels, unless they 
were accompanied by a corresponding decrease in the radial transport timescale, r. Such 



an mverse re 



ationship is expected if t 



interchange (ISouthwood and Kivelson 



le radial transport of plasma is driven by flux tube 



1989 



Brown and Bouchez 



1997 



Pontius et al 



19981). 
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Since, the exact functional form of the inverse relationship is dependent on the details of 
radial convective transport, which remain poorly understood, a variety of inverse relations 
were tested and a rela tion of r oc was adopted. In this work, we extend the model of 



Delamere et al. 



(120041 ) to include an azimuthal dimension. 



1.2. Azimuthal variability during the Cassini era 



Superimposed on the long-term temporal variation in torus composition, UVIS 
observed a persistent azimuthal, i.e., longitudinal, asymmetry in plasma compo sition, 



20061 ). This 



electron temperature, and equatorial electron column density (iSteffl et al. 
nearly-sinusoidal azimuthal variation can be clearly seen in Fig. [TJ The azimuthal variations 
of S II, S III, and electron column density mixing ratios are all approximately in phase with 
each other and are approximately 180° out of phase with the variations of the mixing ratios 
of S IV and O II and the torus equatorial electron temperature. 



EDITOR: PLACE FIGURE □ HERE. 



Over short timescales (<50 hours), the observed azimuthal variability in the lo plasma 
torus is well described by a simple sinusoidal curve with a period equal to the 9.925-hour 
System III rotation period. However, the observed phase of the sinusoidal variation slowly 
drifts to greater System HI longitudes, at a rate of 12.5°/day. This effect is clearly seen in 
the upper panel of Fig. [2l The rate of the phase increase implies that the compositional 
variations observed by UVIS rotate Jupiter with a period of 10.07 hours, 1.5% longer than 



the Sy stem HI rotation pe r iod, y e t 1.3% s 



period 



Sandel and Desslen (119881 ): 



Brown 



using Lomb-Scargle periodograms (ILomb 



lorter than the previously observed "System IV" 



1976 



19951) . Careful ana l ysis of the UVIS observati ons 



Scargle 



1982 



Horne and Baliunas 



19861 ) 



confirmed strong torus periodicity with a period of 10.07 hours (see Fig [7]) and a secondary 
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periodicity with a period close to the System III period of 9.925 hours. 



EDITOR: PLACE FIGURE [2] HERE. 



Numerous ground-based observations have established that plasma in the lo torus la gs 



rigid co-rotation with 



1994a 



Thomas et al. 



Jupiter's magnetic field (iBrown 



198^ 



Roesler et al 



1984 



Brown 



200ll ). It is therefore tempting to think of the 10.07-hour periodicity 



as direct evidence of subcorotating plasma. However, spectroscopic measurements of 



rotation speed of torus plasma have shown that t 



a function of radial distance from Jupiter (IBrown 



le co-rotational lag varies s ignificantly as 



1994a 



Thomas et al. 



200ll ). whereas the 



phase drift seen in the UVIS d ata remains c oherent over a wide range of radial distances. A 



similar argument was used by iBrownl (1l995l ) to rule out plasma subcorotation as the cause 
of the 10.21-hour "System IV" periodicity seen in ground-based optical and Voyager 1 
UVS observations of the lo plasma torus. Instead, the 10.07-hour periodicity in the UVIS 
data, which is phenomenologically similar to the 10.21-hour "System IV" periodicities. 



appears to be t 



le result o: 



( iBrown 



1994b 



Steffi 



2005 



a compositional w ave propagating azimuthally through the torus 



Steffl et al. 



20061 ). Given the apparent simil arities between the 



0.0 7- hour periodicity in the UVIS data, t 



le 10. 2 2 4- hou r periodicity of 



( 119881 ). and the 10.214-hour periodicity of 



Sandel and Dessler 



Brownl (119951 ). we will subsequently refer to all 



three phenomena as "System IV". 

In addition to the phase drift, the relative amplitudes of the azimuthal variations in 
S II and S IV mixing ratios vary, in a roughly cyclical manner, between 5-25%, as shown in 
the lower panel Fig. [2l The time between the observed peaks of the S II and S IV azimuthal 
variations is ~29 days, suggestively close to the 28.8-day beat period between the observed 
10.07-hour Cassini epoch System IV period the 9.925-hour System HI rotation period. 
Thus, the amplitude of the azimuthal variation of these two ion species appears to be 
modulated by the its location, relative to System HI longitude. The amplitude is greatest 
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when the peak S II mixing ratio is located near A/77=210±15° and smallest, i.e., most 
azimuthally uniform, when the peak S II mixing ratio is located near A///=30±15°. This 
effect is clearly seen in S II and S IV. However, the amplitude of the azimuthal variation 
in the two primary ion species, O II and S III, was relatively constant during the UVIS 
observing period (in the range of 2-5%). 



2. Azimuthal Models 

The one-box model of DB03 has five input parameters that can be adjusted to match 
the conditions observed in the lo torus: the rate of neutral atoms supplied to the torus 
(Snfl), the ratio of oxygen to sulfur in the supplied neutrals (O/S), the fraction of the 
total electron population heated to super-thermal levels (//i,o), the temperature of the hot 
electrons (Te^hot), and the radial transport time scale (r). The sensitivity of the model to 
these parameters is discussed in detail by DB03. 

A successful azimuthal model must produce variations in torus composition that i) are 
single-peaked and nearly sinusoidal, ii) exhibit the temporal changes seen by Cassini UVIS 
over the 45-day approach phase, iii) drift, relative to System III longitude, at a rate of 
approximately 12.5° day^^, iv) vary in amplitude with a period of 28.8 days, and v) have 
relative phases and amplitudes consistent with those obse rved by UVIS. To me et these 
requirements, we developed 1-d variants of the DB03 and 



Delamere et al 



mQ4 \ models. 



which are discussed below, in order of increasing complexity. 



2.1. Basic Azimuthal Model 



In the basic azimuthal model, which serves as the basis for all subsequent models, we 
extended the one-box model of DB03 to include 24 azimuthal bins, each corresponding to 
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a 15° segment in System III longitude, located at a radial distance of 6 Rj. The density 
and temperature of torus plasma and neutral species are calculated for each azimuthal bin, 
as are reaction rates between species. The model plasma lags co-rotation with the System 
III coordinate system by an amount, Af , which can be a function of System III longitude. 
Several models that incorporated an azimuthally variable co-rotational lag were tested, but 
none could satisfactorily match the UVIS observations. Therefore, we hold Av constant 
with location in the torus. Spectroscopic observations of the lo torus have shown that torus 
plasma at a radial distance of 6 Rj l ags co-rotation with the magnetic f ield (i.e., the System 



III coordinate system) by 3-4 km/s (IBrown 



1994a 



Thomas et al. 



20011); we therefore adopt 



a value of Af=3.5 km/s, corresponding to a rotation period of 10.41 hours.. 

lo's extended neutral clouds are on Keplerian orbits of Jupiter. At QRj, they have a 
velocity of ~57 km s~^ relative to the System III coordinate frame. Thus, neutrals will 
cross a 15° azimuthal bin in ~2000 s. This is much faster than the characteristic timescales 
for torus chemistry (ionization, charge exchange, etc.), and it places an upper limit on the 
length of the model time step. In the model, when a neutral atom becomes ionized, it is 
given the plasma subcorotation velocity Av and pickup energy of 380 eV for sulfur ions and 
190 eV for oxygen ions. 

The creation of lo's extended neutral clouds an d their three-dimensional spatial 



distribution are wel 



2003 



Thomas et al. 



beyo nd the scope of this model (ISmyth and Marconi 



2003 



Burger 



20041 ). Instead, the model starts with an azimuthally uniform neutral 



cloud. New oxygen and sulfur atoms are added in an azimuthally uniform manner at a 
rate controlled by the source parameters in Eq. HJ Snfl, a„, t„, and (T„. The rate at which 
neutral atoms are lost to the torus, however, will generally vary as a function of azimuthal 
position, resulting in neutral density variations of up to 20%. We assume the neutrals 
are confined to Jupiter's rotational equator with a scale height of 0.5 Rj, consistent with 
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Burgerl (120031 ): our models results are insensitive to changes in the neutral scale height. 



The transport of mass and energy between azimuthal bins is handled via a two-step 



Lax-Wendroff scheme (IPress et al. 



19921 ) implemented after the densities and temperatures 
for have been updated. Following the azimuthal transport, sinusoidal curves are fit to the 
azimuthal variations of the model, in a procedure equivalent to that used in the analysis of 
the UVIS data (cf. Steffi et al., 2004). 

Torus plasma experiences a significant centrifugal force due to the rapid rotation of 
Jupiter and finds an equilibrium about the position on a given magnetic field line that is 
most distant from Jupiter's rotation axis (see Bagenal, 1994). The locus of these points 
forms the centrifugal equato r, which is locat e d 1/3 of the way betwe en the magnetic equator 



and the rotational equator (IHill et al. 



1974 



Cummings et al. 



19801 ). The offset between 



the centrifugal equator and Jupiter's rotational equator varies as a function of System III 
longitude, ranging from Rj at A///=20° and A///=200° to 0.67 Rj at Am=110° and 
A///=290°) at a radial distance of 6 Rj. 

Plasma pressure forces cause the torus plasma to spread out from the centrifugal 
equator along magnetic field lines wit h a scale height determined by the mass and 



temperature of the ions (IBagenal 



range between 1-2 Rj (ISteffi 



994j ). Since typical scale heights of ions in the lo torus 
20051 ). the number density of torus ions at the rotational 
equator can vary by up to ~ 40%. Thus, timescales for torus reactions are affected by the 
latitudinal density distribution, with ion/neutral interactions affected more strongly than 
ion/ion interactions. 

Although our model includes only one spatial dimension (azimuthal), it includes the 
effects of the latitudinal distr i butio n of torus plasma via the method of latitudinal averaging 



described by 



Delamere et al 



(I2OO5I ). This method approximates the latitudinal distribution 



of ions (neutrals) as a Gaussian centered on the centrifugal (rotational) equator. The 
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temperature distribution of each species is assumed to be Maxwellian, so that temperature 
remains constant with latitude. Reaction rates for each azimuthal bin are determined by 
calculating the average reaction rate weighted by the latitudinal distribution of the two 
react ants. 

The DB03 model includes only the five major ion species in th e torus: S II, S III, S IV, 



O II, and O III. Several oth er min or species such as CI II and CI II 



2000 



Feldman et al 



20011 ). C III flFeldman et al. 



2004h . and S V f lSteffl et al. 



(iKuppers and Schneider 



2004br ) have 



been detected in the torus. Additionally, the presence of Na and K ions is inferred, given 
the observations of neutral Na and K near lo (see review by Thomas et al., 2004). Finally, 
the torus plasma will also contain protons. Alt hough early work estimated the total flux 



tube content consisted of 10-15 % protons, e.g., ( 



this value at only a few percent (ICrary et al 



To 



1996 



car et al. 



Wang et al 



19821). mo re recent work limit s 



1998a 



b 



Zarka et al. 



200ll ). 



Although we do not include any reactions involving these minor species, they are included 
in our model's calculation of charge neutrality by assuming that they compose 10% of the 
total charge, such that: 



0.9Nr = 



2S' 



3S^ 



0+ + 20' 



(2) 



A similar equation for quasi-neutrality was used by ISteffl et al.l ( 120061 ) in the analysis of 
the Cassini UVIS spectra. Model results are generally insensitive to the assumed charge 
fraction of these minor ion species, over the range of 0-20%. 

To accommodate the variations to the basic azimuthal model described below, the 
fraction of hot electrons was allowed to vary as a function of both time and location 
according to the following equation: 
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(3) 



X (1 + ah,\jy COs(A/77 - 4>h,\iv - ^t) + (^h,\ni COs(A/// - (phMn)) 



However, in the basic model, ah, ahXjy, and ahx,,, are set to zero. 



2.2. Time- Variable Model 



To reproduce the temporal 
epoch, we follow the method of 
the neutral source rate: 



changes in toru s comp osition observed during the Cassini 



Delamere et al. 



(120041 ) and include a Gaussian increase in 



Sn{t) = Sn,0 (^l + 



(4) 



where Snfl is the baseline value of the neutral source rate. Like iDelamere et al.l (120041 ) . the 
date of the peak neutral source, t„, was held fixed at day 249 to match the center of the 



incre ase in logenic dust flux observed by the Galileo Dust Detector System (iKriiger et al. 
20031 1. 



Latitudinal averaging of reaction rates produces an increase in the characteristic 
timescales for torus reactions. As a result, although the general trend of decreasing S II 
and increasing S IV seen in the UVIS data (c.f. Fig. H]) could be reproduced, the rapidity 
of this change could not. Based on their parameter sensitivity studies, DB03 noted that 
the hot electron fraction {fh) is the only model parameter capable of modifying the torus 
composition on such rapid timescales (days). Therefore, a broad Gaussian perturbation to 
the hot electron fraction was added to the model by allowing ah to be non-zero, th-, and ah 
were also allowed to vary to match the observed compositional change. 
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Adding a Gaussian increase in the hot electron fraction to match the UVIS-observed 
composition is admittedly ad hoc. However, this increase can be loosely justified by the 
following argument. As the the density of the neutral clouds increases, there is increased 
mass loading in the torus. The resulting increase in fiux tube content, increases the 
efficiency of outward radial transport of plasma, parameterized in our model as r. This, in 
turn, will drive field-aligned currents that could produce additional hot electrons. 



2.3. Subcorotating Hot Electron Model 



The azimuthal variation in torus composition observed by Cassini UVIS has a rotation 
period of 10.07 hours, slightly longer than the 9.925-hour System III rotation period. 



The torus plasma, however, has an even longer ro tation period (roughly 10.41 



6 Rj), and is a strong function of radial distance (IBrown 



lours at 



1994a 



Thomas et al. 



200l|), so 



the subcorotation of the torus plasma can not be directly responsible for the period of 
the azimuthal variations. Instead, noting that the timescale for hot electrons to couple 
energetically with the thermal electron population is of order tens of minutes compared to 
timescales of several days or more for all other torus processes (cf. Table E]), we introduce a 
subcorotating variation in the fraction of hot electrons by allowing ah^Xjv and (ph,Xiv Eq. [3] 
to be non-zero. The angular velocity of the hot electron variation relative to System III, 
u, is set to 12.5°/day, corresponding to a rotation period of 10.07 hours. We subsequently 
refer to this model as the "Subcorotating Hot Electron Model" . 



2.4. Dual Hot Electron Model 



In addition to having a period of 10.07 hours, the amplitude of the azimuthal variation 
in torus composition appears to change with time in a roughly periodic way. The period 
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of the amplitude variation (28.8 days) is identical to the beat period produced by the 
interference of the observed 10.07-hour period and the System III period of 9.925 hours. 
A similar variation in the UV brightness of the torus was observed by the Voyager 2 UVS 



( ISandel and Dessler 



19881 ): the 14.1-day period of this brightness variation matches the beat 



period between the Voyaqer-eia. System IV period of 10.22 hours and the System III period. 



Further analysis by lYang et al. 



(119911 ) concluded that the System IV periodicity in the 
Voyager daia, was indeed independent of System III and that the observed 14.1-day period 
was likely the result of the interference of phenomena with the System IV and System III 
periods. 

To test whether a similar beat frequency modulation of torus composition could be 
produced by super-thermal electrons, we added a second hot electron variation that remains 
fixed in System III coordinates by allowing (yh,\iii to be non-zero. This model, with all 
sixteen input parameters non-zero, will be referred to as the "Dual Hot Electron Model" . 



Model Results 



3.1. Basic Azimuthal Model 



The five parameters of the basic azimuthal model {Snfi, 0/S, fhfl, Te^hot, and r) were 
adjusted to match the azimuthally-aver aged torus compo sition derived from the Cassini 



UVIS observations of 2001 January 14 (jSteffl et al. 



2004bl ). The values of these parameters 



are shown in Table [H The final equilibrium state produced using these parameters is used 
as the initial condition of the lo torus for the subsequent three models. 



EDITOR: PLACE TABLE □ HERE. 
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Characteristic timescales for torus loss processes are given in Table [2l which illustrates 
the importance of charge exchange reactions (and resonant charge exchange reactions 
in particular) in the lo plasma torus. For example, while the primary loss mechanism 
of neutral sulfur from the extended clouds is electron impact ionization by the thermal 
(~5 eV) electron population, neutral oxygen is lost primarily through resonant charge 
exchange with O II. Likewise, the dominant loss process of S II is not ionization (by either 
the thermal or hot electron population), but rather the resonant charge exchange reaction 
+ S^^ S^^ + . Although resonant charge exchange reactions do not change the 
number densities of torus ions, they do redistribute energy between ion species. 

EDITOR: PLACE TABLE M HERE. 

The offset between the centrifugal and rotational equators produces a slight azimuthal 
variation in torus composition, as can be seen in Fig. [31 Where the two equator planes 
intersect (at 110° and 290°), S II shows a 1% increase, relative to the average mixing ratio, 
due primarily to the increased rate electron impact ionization of S I. Approximately 30° 
downstream (torus plasma moves in the direction of increasing System III longitude), S IV 
exhibits a 1% decrease due primarily to the increased rate of the charge exchange reaction 
5^+ + O ^ 5^+ + O^. However, this azimuthal variation is double-peaked and clearly much 
smaller in amplitude than the azimuthal variations observed by Cassini UVIS. 

EDITOR: PLACE FIGURE [3] HERE. 

3.2. Time- Variable Azimuthal Model 

After including Gaussian perturbations to both the neutral source rate and the hot 
electron fraction, the time-variable azimuthal model reproduces the observed temporal 
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behavior of the sulfur ion mixing ratios, as seen in Fig. HI The parameter values used 
to produce this figure are listed in Table [TJ However, the model fails to reproduce the 
temporal behavior of O II. This discrepancy could arise from an error (or errors) in the 
rate coefficients for reactions involving O II. Alternatively, the ratio of oxygen to sulfur 
atoms sup plied to the extendec 



possibility, iSpencer et al. 



neutral clouds may not be constant. Regarding this latter 
(120001 ) report the discovery of gaseous S2 in Jo's Pele plume at 
the level of SO2/S2 = 3-12. If the putative volcanic event responsible for the increase in 
the torus neutral source (the Tvashtar eruption of 2000) was sufficiently rich in S2 (or SO), 
the ratio of oxygen to sulfur atoms supplied to the neutral clouds could have temporarily 
decreased. As the neutral source rate returned to pre-event levels, the 0/S ratio of the 
neutrals would rise, producing a gradual increase in the mixing ratio of O II. 

EDITOR: PLACE FIGURE H HERE. 



3.3. Subcorotating Hot Electron Model 

As expected, the subcorotating sinusoidal variation in the hot electron fraction 
produced single-peaked azimuthal variations in torus composition with a rotation period of 
10.07 hours. The top panel of Fig [5] shows the model azimuthal variation drifts in System 
III longitude by 12.5° per day and roughly matches the observed phase increase. The 
best-fit parameter values of the subcorotating hot electron model are listed in Table [H 

Including an azimuthal variation in the hot electron fraction did not change the average 
composition of the torus. This behavior, due to the symmetric nature of the azimuthal 
perturbation (a sine wave), greatly simplifies the fitting procedure. Since the temporal 
variation model parameters (a„, t„, cr„, ah, t^, and ct/^) are decoupled from the azimuthal 
variation model parameters (a/i,Ajv 4>h,Xjv)^ once the best-fit temporal parameters have 
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been determined only two additional parameters need to be varied to match the phase 
increase of the azimuthal variation in composition. 

The interaction of the azimuthal variation in hot electron fraction with the neutral 
source increase results in a dramatic increase (from 6% to 22%) in the relative amplitude 
of the S IV azimuthal variation, centered around DOY 249. This increase is largely caused 
by the efficient removal of S IV ions via the charge exchange reaction S^^ + O ^ 5*^+ + 
at longitudes where the centrifugal and rotational equators intersect. During the Cassini 
UVIS approach observations (DOY 275-320), the amplitude of the model S IV azimuthal 
variation monotonically decreases to its pre-event level. Conversely, interaction with the 
increased neutral source results in a slight decrease (from 12% to 10%) in the amplitude 
of the model S II azimuthal variation. During the Cassini UVIS observation period, the 
S II variation exhibits a slight maximum near day 295. The amplitude behavior of both 
the S IV and S II variations produced by the subcorotating hot electron model is in stark 
contrast to the UVIS observations, as seen in the bottom panel of Fig. [51 

EDITOR: PLACE FIGURE [5] HERE. 

3.4. Dual Hot Electron Model 

Like the subcorotating hot electon model, the dual hot electron model can produce 
azimuthal variations in torus composition that drift to higher System III longitudes, 
but only if ah,\iv ^ at^Xjjj. If this condition is not satisfied, the azimuthal variation in 
composition becomes fixed in System III coordinates. In contrast to the subcorotating 
hot electron model, the rate at which the azimuthal variation drifts is not constant. This 
effect can be seen in the top panel of Fig. [61 For most of the time, the model azimuthal 
variation lags co-rotation with a period of 10.02 hours, corresponding to an angular velocity 
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in the System III frame of u = 8.3°/day. However, when the two variations are roughly 
180° out of phase, the rotation period of the compositional variation increases to 10.23 
hours {uj = 26.1°/day). The UVIS data show a similar trend, though the difference in 
rotation period is less dramatic, changing from 10.02 hours {u = 8.8°/day) to 10.11 hours 
{uJ = 16.2°/day). The changing rotation period of the compositional variation is in no way 
related to the speed of the torus plasma, but rather is a wave phenomenon produced by the 
interaction of the two hot electron variations with the torus plasma. 

EDITOR: PLACE FIGURE [6] HERE. 

Although the instantaneous rotation period of the azimuthal variation varies over the 
28.8-day beat period, Lomb-Scargle periodogram analysis shows strong periodicity at 10.07 
hours, with secondary periodicity at the System HI period of 9.925 hours, just like the 
UVIS data (see Fig. [7j). For a stable neutral source (a„ = ah = 0) and the Cassini epoch 
nominal torus composition, the amplitude of the S II variation ranges from 5% to 17% over 
the 28.8-day beat period. For S IV, these values are 3% and 10%. 

EDITOR: PLACE FIGURE [7] HERE. 

As with the subcorotating hot electron model, the interaction of the neutral source 
increase with the two azimuthal variations in hot electron fraction results in a dramatic 
increase in the amplitude of the variation of S IV (7% variation at minimum amplitude 
to 32% variation at maximum amplitude), while producing only a minimal effect on the 
amplitude of the S II variation. As the neutral source rate gradually returns to its pre-event 
level, the amplitude of the S IV variation also decreases, reaching its pre-event amplitudes 
after approximately 70 days. Similar behavior can be seen in the amplitude of the S IV 
variation seen by UVIS (see bottom panel of Fig. [6]). 
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Figure [8] shows the azimuthal variation in torus composition produced by the dual 
hot electron variation model during epochs of maximum (Day 279) and minimum (Day 
293) amplitude compared to the azimuthal variation observed by Cassini UVIS. The model 
parameters used to produce the results shown in Figs. [7H8] are given in Table [H 

EDITOR: PLACE FIGURE E] HERE. 
4. Discussion 

Given the simplifying assumptions, the match between the output of the dual hot 
electron variation model and the Cassini UVIS observations of the lo torus is remarkable. 
Like the observed variations, the model variations are single-peaked and have an average 
rotation period of 10.07 hours. The model variations show the same relative phase and 
amplitude as the observed variations, exhibit a 28.8-day beat period, and, with the 
exception of O II, match the azimuthally-averaged torus composition as a function of time. 

The large number of model input parameters (16), makes it difficult to assess the 
uniqueness of our model solution, as a complete search through the 16-dimensional 
parameter space is computationally prohibitive. However, despite considerable effort, we 
found no other region in the parameter space of our model that produced results that match 
the UVIS observations. Regarding the uncertainty in parameter values, models run with 
deviations of up to 10% from the parameter values given in Table [1] generally yield results 
that are similar to those shown in Figs. [6] and [H whereas models with parameter deviations 
greater than ~10% produce a notably poorer match to the UVIS data. 

The value of the model parameter (pjy provides no real insight into the nature of the 
torus; rather, this parameter merely represents the phase of the subcorotating hot electron 
variation at an arbitrary time, t = (in this case, chosen to be 2000-01-01). In contrast. 
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the requirement that = 290° for the model to match the UVIS data is significant. 
Since Jupiter's dipole magnetic field is tilted toward a System III longitude of 200°, the 
rotational and centrifugal equator planes intersect at the System III longitudes of 110° and 
290°. In the absence of any temporal variations, this results in a 20% increase in the pickup 
energy supplied to the lo torus at these longitudes. By itself, this has very little effect on 
the torus composition, as shown by Fig [31 If, as we propose, the super-thermal electrons in 
the lo torus are primarily produced via field-aligned currents, the increase in mass loading 
at these longitudes should drive additional field-aligned currents resulting in an increase of 
hot electrons at these longitudes. However, increased mass loading (and thus an increased 
hot electron fraction) should occur at both A/// = 290° and Xui = 110°. Models using 
a double-peaked System Ill-fixed hot electron variation were not able to reproduce the 
amplitude variations seen in the UVIS data and resembled the results of the Subcorotating 
Hot Electron model shown in Fig. [51 Our best model requires a hot electron maximum 
at A/// = 290°, as expected, but a hot electron minimum at A/// = 110°. What causes 
the symmetry between Xm = 110° and A/// = 290° to be broken? One possibility is that 
higher-order components of Jupiter's magnetic affect the conductivity of the ionosphere in 
such a way that the field-aligned currents that accelerate electrons in the torus are strongly 
favored at Am = 290° over Xm = 110°. 

While the System Ill-fixed variation in hot electrons likely has its origin in the 
interaction between Jupiter's magnetic field and the Jovian ionospher e, the source o f the 
System IV periodic i ty in the lo torus has been enigmatic. Previously, 



Desslen (Il985l ) and 



Sandel and Desslerl (Il988l ) have proposed a secondary, high-latitude component of Jupiter's 
magnetic field that lags co-rotation by a few percent could be responsible for producing 
the System IV periodicity. While neither the UVIS observations nor our efforts to model 
them can rule out the existence of such a high-latitude component, this hypothesis is 
problematic. Presumably, such a magnetic field component would affect the Jovian aurora 
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as well as the lo plasma torus. However, the morphology of the Jovian aurora is strongly 
fixed in System III longitude, and although the temporal coverage of the Jovian aurora 
has been somewhat sporadic, there have been no aur oral phenomena rep orted at either 



the 10.2 or 10.07-hour periods typified by System IV (j Clarke et al. 



20041 ). Furthermore, it 



is difficult to reconcile a secondary magnetic field component wi th the changing period o f 



System IV phenome na: 10.2 2 4 hou rs during the Voyager epoch (jSandel and Desslei 



10.214 hours in 1992 



Brownl (Il995l ) and 10.07 hours during the Cassini epoch (ISteffl et al 



1988h. 



20061 ). We therefore consider the existence of a subcorotating high- latitude magnetic field 
component improbable, though we are unable to offer a satisfactory alternative. 

As discussed above, the rotation speed of torus plasma is a strong function of radial 
distance from Jupiter. Both the Subcorotating Hot Electron Model and the Dual Electron 
Model use a corotational lag of 3.5 km/s. However, both models are insensitive to changes 
in the amount of corotational lag, over the range observed in the torus (0-4 km/s). In 
particular, the Dual Hot Electron Model can reproduce the temporal and azimuthal 
variations in composition observed by UVIS, regardless of the radial distance at which it 



i s run , consistent with the radial uniformity of the System IV period observed by 



Brown 



5. Conclusions 

During the Cassini spacecraft's fiyby of Jupiter, the UVIS instrument observed 
remarkable temporal and azimuthal variations in the composition of the lo plasma torus. 
The azimuthal variations, which are primarily seen in the ion species S II and S IV, lag 
co-rotation with the magnetic field and are decoupled from the rotation speeds of both the 
torus plasma and logenic neutral clouds. The strength of the azimuthal variation changes in 
a seemingly periodic manner with a period of 28.8 days — the beat period between System 
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III and the observed 10.07-hour rotation period of the azimuthal variations. 

To model the tempo ral and azimuthal chang e s obse rved by UVIS, we have extended the 



torus chemistry model of 



Delamere and Bagenall (120031 ) to include an azimuthal dimension. 



Our preferred model, which includes two independent azimuthal variations in the amount of 
hot electrons in the lo torus, one subcorotating and one fixed in System III, can reproduce 
the UVIS observations remarkably well. The major findings of this paper are summarized 
below. 



1. The months- long change in the average composition of the lo plasma torus can be 
modeled by introducing a factor of 3 increase to the rate of oxygen and sulfur atoms 
supplied to the extended neutral clouds that are the source of the torus plasma 
coupled with a 30% increase in the fraction of hot electro ns in the lo torus. This 



result is similar to that reported by 



Delamere et al 



fl2004h . 



2. An azimuthal variation in the fraction of hot electrons in the lo torus that rotates 
with a period of 10.07 hours can produce subcorotating azimuthal variations in torus 
composition like those observed by Cassini UVIS. 

3. The interference of the subcorotating hot electron variation with a second hot electron 
variation that remains fixed in System III can produce the beat frequency modulation 
in the amplitude of the azimuthal variations also seen by Cassini UVIS 



A. Model Reaction Rate CoefRcients 



We describe here the sources of the various rate coefficients used in our model. 
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A.l. Ionization 



Rate coefficients for electron impact ionization are calculated using the fit formulae 
given by 



Voronov 



recommended data (IBell et al 



9971). These fi t s are based on the University of Belfast group 



1983 



Lennon et al 



19881). 



Since the lifetimes of neutral and ionic oxygen and sulfur against photoionization are 
several orders of magnitude longer than the characteristic timescales for other processe s 



such as electron impact ionization, charge exchange, and recombination (iHiibner et al. 
19921 ). the effects of photoionization can be ignored. 



A. 2. Recombination 



For ions with multiple electrons, recombination rate coefficients are usually divided into 



two s eparate processes: radiative recombination and dielectronic recombination (jOsterbrock 



19891 ). The total recombination rate coefficient is the sum of the radiative and dielectronic 
recombination terms. Tota l recombination rate coefficients for oxygen ion species are 

). Total recombination rates for S^^ 



obtained from 
obtained from 



Nahai 



Naharl ( 



(11999 



S2+ and S2+ ^ S+ are 



19951 ) and the associated erratum iNaharl (119961 ). In general, the 



rate coefficients published by Nahar agree well with previously published results at low 
temperatures. At temperatures typical of the lo torus, however, the Nahar rates can be up 
to an order of magnitude lower than previously published values. Since the work by Nahar 
is the most recent treatment of the recombination rate problem for sulfur and oxygen ions 
and employs a more sophisticated technique than previous studies, these rates are used in 
the model. For the recombination of S+ — > S, the radiative recombination rate coefficient of 



ShuU and van Steenberd ( 



is obtained from 



19821) i s used , while the dielectronic recombination rate coefficient 



Mazzotta et al. 



(11998h . 
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The recombination rate coefficient s used in this work differ from those used in previou s 



versio ns of the torus chemistry model (jPelamere and Bagenal 



2003 



Delamere et al. 



2004 



20051 ). However, since recombination reactions are generally much slower than other 
processes that occur in the torus (cf. Table [2]) these changes do not significantly affect our 
conclusions. 



A. 3. Charge Exchange 



Charge exchange reactions play an important role in the torus chemistry. Seventeen 
charge exchange reactions between atomic and ionic species of sulfur and oxygen, listed in 
Tabl e 1 of DB03, are includ e d in t he model. All charge exchange reaction rates are taken 



from 



McGrath and Johnson] fll989h . 



A. 4. Radiation 



The radiative rate coefficients of the model are obtained from the CHIANTI atomic 



physics database version 4.2 (jPere et al. 



1997 



Young et al 



20031). 
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Table 1. Best-fit Azimuthal Model Parameters 



Parameter Description 


Symbol 


Basic 


Time Variable 


Subcorotating Hot 


Dual Hot 








iviouci 






Neutral source rate at the rotational equator (cm~^ 


Sn,0 


9.9x10-* 


9.9x10"* 


9.9x10-* 


9.9x10"* 


O to S ratio of neutral source 


O/S 


1.55 


1.55 


1.55 


1.55 


Fraction of hot electrons 


fh,0 


0.00235 


0.00235 


0.00235 


0.00235 


Temperature of hot electrons (eV) 




55 


55 


55 


55 


Radial transport timescale (days) 


TO 


62 


62 


62 


62 


Amplitude of neutral source increase 


CXn 


0.0 


2.4 


2.4 


2.4 


Date of neutral source increase (DOY 2000) 


tn 




249 


249 


249 


Gaussian width of neutral source increase (days) 






30 


30 


30 


Amplitude of hot e~ fraction increase 




0.00 


0.30 


0.30 


0.30 


Date of hot e~ fraction increase (DOY 2000) 


th 




279 


279 


279 


Gaussian width of hot e~ fraction increase (days) 






60 


60 


60 


Amplitude of System IV hot e~ variation 


0'h,Xiv 


0.0 


0.0 


0.43 


0.43 


Phase of System IV hot e~ variation (°) 


<t'h,Xlv 






60 


60 


Angular velocity between Systems III and IV (°/day) 








12.5 


12.5 


Amplitude of System III hot e~ variation 




0.0 


0.0 


0.0 


0.30 


Phase of System III hot e~ variation (° ) 


<t>h,XijI 








290 
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Table 2. Characteristic Timescale of Torus Loss Processes 



Loss Mechanism 


S I 


S II 


S III 


S IV 


O I 


O II 


O III 


Thernial e impaiCt ionizaition 


U.o 


1 A n 

lO.U 




1 riztnn 






( u / uu 


Hot e impact ionization 


lo.y 




1 OS 


OOCS 


"4:0.0 


1DC5 




R.6Combinat ion 




1 /1 1 n 




1 O'i 




^tuou 


1 QQfl 

iooU 


G ' -ho'' — > o ' ' -ho' 




o.u 












o "h o — * o ~r o 


o.u 














Q _i_ + + - C+ 1 Q + 

D -r o — o ~r o 


















4.0 




240 












14.0 






142 








o + o+ o+ + o* 










2.6 


43.3 




o + 0++ o+ + o+ 










627 




1070 


O + 0++ ^ 0++ + o* 










60.4 




104 


o + s+ ^ o+ + s* 




8510 






1990 






S + 0+ ^ s+ + o* 


10.8 










734 




S + 0++ ^ S+ + 0+ 


13.9 












95.7 


S + 0++ S++ + 0+ + e- 


20.1 












138 


O + S++ -* 0+ + S+ 






205 




13.8 






0++ + S+ -» 0+ + S++ 




262 










105 


+ S+++ ^ 0+ + S++ 








24.4 


9.6 






0++ + S++ ^ 0+ + S+++ 






376 








43.4 


S+++ + S+ ^ S++ + S++ 




585 




346 








Radial transport 




62.0 


62.0 


62.0 




62.0 


62.0 


Total of all loss processes 


0.5 


2.2 


7.3 


12.8 


1.3 


20.9 


12.5 



Note. — Characteristic timescales given are for the "nominal" Cassini epoch torus com- 
position and have been aaimuthally averaged. All timescales have units of days. 
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Fig. [T] Caption 

Relative ion mixing ratios, electron temperature, and electron column density for a typical 
3-day period obtained from the dusk ansa. Values are normalized to the average value over 
the 3-day period. The best-fit sinusoids for this period are overplotted. Note the strong 
anti-correlation of S II with S IV and equa torial electron tem perature with equatorial 



electron column density. Figure taken from iSteffl et al. 



(120061). 



Fig. H Caption 

Azimuthal variations in the lo plasma torus as observed by Cassini UVIS. The top panel 
shows the location (in System III coordinates) of the peak in mixing ratio of the primary 
ion species in the lo torus as a function of time. All four ion species show a roughly linear 
trend of increasing phase with time. The bottom panel shows the relative amplitude (as a 
percentage) of the azimuthal variations as a function of time. The relative amplitudes of 



O II and S III remain around the few percent level, while the con iparatively 



ion species S II and S IV vary between 4-25%. Figure taken from iSteffl et al. 



ess abundant 



(120061). 



Fig. M Caption 

Relative variation in torus composition produced by the basic azimuthal model. The 
small azimuthal variations are caused by the offset between the centrifugal and rotational 
equators and are therefore double-peaked. 



Fig. m Caption 

Azimuthally-averaged composition of the lo torus as observed by Cassini UVIS and 
reproduced by the time-variable azimuthal model. Observed mixing ratios derived 
from UVIS spectra are shown with plot symbols and black connecting lines. Uncer- 
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tainties in the UVIS-derived mixing ratios are approximately 10%, as shown by the 
error bars. Averaged mixing ratios produced by the time-variable azimuthal model are 
shown with thick solid lines. Although the model results shown here were produced by 
the time-variable azimuthal model, virtually identical results can be produced by the 
subcorotating hot electron model (Section [23l) and the dual hot electron model (Section [2^ . 

Fig. [5] Caption 

Comparison of Cassini UVIS data with output from the subcorotating hot electron model. 
UVIS observations of both dawn and dusk ansae have been averaged together. The top 
panel shows the azimuthal location (in System III coordinates) of the peak mixing ratios 
of S II and S IV. The bottom panel shows the amplitude of S II and S IV azimuthal variations. 

Fig. [6] Caption 

Comparison of Cassini UVIS data with output from the dual hot electron model. This 
model features the superposition of two azimuthal variations in hot electron fraction, 
one with a rotation period of 10.07 hours and the other with the System III rotation 
period of 9.925 hours. Unlike the azimuthal variations produced by the subcorotating hot 
electron model (Fig [5]) the phase of the azimuthal variation increases more rapidly when 
the amplitude is near its minimum value (near DOY 293). The bottom panel shows the 
amplitude of azimuthal variations of S II and S IV. The interference of the two hot electron 
variations creates a beat period of 28.8 days. 

Fig. [7] Caption 

Lomb-Scargle periodograms of the mixing ratio of S II as derived from Cassini UVIS data 
(left) and the dual hot electron model (right). Model mixing ratios have been sampled 
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at the same time and spatial location as the UVIS observations. Both data and model 
periodograms show a sharp peak at a frequency of 0.0993 h"^, corresponding to the 
10.07-hour "System IV" period observed by Cassini UVIS and a secondary peak at the 
System III rotational frequency. Sidebands of the two peaks can be seen near 0.05 h~^ and 
0.20 h~"^. Small spurious peaks due to the sampling interval of the UVIS data are also present. 

Fig. il Caption 

Azimuthal variation of ion mixing ratios in the lo plasma torus during two 2-day periods. 
Plotting symbols represent ion mixing ratios derived from Cassini UVIS data: diamonds 
from the dawn ansa and pluses from the dusk ansa. The solid lines are output from the 
dual hot electron model. 
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Fig. 1.— Steffi et al. 

Relative ion mixing ratios, electron temperature, and electron column density for a typical 
3-day period obtained from the dusk ansa. Values are normalized to the average value over 
the 3-day period. The best-fit sinusoids for this period are overplotted. Note the strong anti- 
correlation of S II with S IV and equatorial electron temperature with equatorial electron 



column density. Figure taken from 



Steffi et al. 



(12006h . 
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Fig. 2.— Steffi et al. 

Azimuthal variations in the lo plasma torus as observed by Cassini UVIS. The top panel 
shows the location (in System III coordinates) of the peak in mixing ratio of the primary 
ion species in the lo torus as a function of time. All four ion species show a roughly linear 
trend of increasing phase with time. The bottom panel shows the relative amplitude (as a 
percentage) of the azimuthal variations as a function of time. The relative amplitudes of 
O II and S III remain around the few percent level, while th e comparatively le ss abundant 



ion species S II and S IV vary between 4-25%. Adapted from 



Steffi et al 



mm . 
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Fig. 3.— Steffi et al. 

Relative variation in torus composition produced by the basic azimuthal model. The small 
azimuthal variations are caused by the offset between the centrifugal and rotational equators 
and are therefore double-peaked. 
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Fig. 4.— Steffi et al. 

Azimuthally-averaged composition of the lo torus as observed by Cassini UVIS and repro- 
duced by tlie time-variable azimutfial model. Observed mixing ratios derived from UVIS 
spectra are shown with plot symbols and the thin connecting lines. Uncertainties in the 
UVIS-derived mixing ratios are approximately 10%, as shown by the error bars. Averaged 
mixing ratios produced by the time-variable azimuthal model are shown with thick solid 
lines. Although the model results shown here were produced by the time-variable azimuthal 
model, virtually identical results can be produced by the subcorotating hot electron model 
(Section 12.31) and the dual hot electron model (Section l2.4p . 
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Fig. 5.— Steffi et al. 

Comparison of Cassini UVIS data with output from the subcorotating hot electron model. 
UVIS observations of both dawn and dusk ansae have been averaged together. The top panel 
shows the azimuthal location (in System III coordinates) of the peak mixing ratios of S II 
and S IV. The bottom panel shows the amplitude of S II and S IV azimuthal variations. 
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Fig. 6.— Steffi et al. 

Comparison of Cassini UVIS data with output from the dual hot electron model. This model 
features the superposition of two azimuthal variations in hot electron fraction, one with a 
rotation period of 10.07 hours and the other with the System III rotation period of 9.925 
hours. Unlike the azimuthal variations produced by the subcorotating hot electron model 
(Fig [5]) the phase of the azimuthal variation increases more rapidly when the amplitude 
is near its minimum value (near DOY 293). The bottom panel shows the amplitude of 
azimuthal variations of S II and S IV. The interference of the two hot electron variations 
creates a beat period of 28.8 days. 
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Fig. 7.— Steffi et al. 

Lomb-Scargle periodograms of tlie mixing ratio of S II as derived from Cassini UVIS data 
(left) and the dual hot electron model (right). Model mixing ratios have been sampled at the 
same time and spatial location as the UVIS observations. Both data and model periodograms 
show a sharp peak at a frequency of 0.0993 h~^, corresponding to the 10.07-hour "System 
IV" period observed by Cassini UVIS and a secondary peak at the System III rotational 
frequency. Sidebands of the two peaks can be seen near 0.05 h~^ and 0.20 h~^. Small 
spurious peaks due to the sampling interval of the UVIS data are also present. 
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Fig. 8.— Steffi et al. 

Azimuthal variation of ion mixing ratios in tlie lo plasma torus during two 2-day periods. 
Plotting symbols represent ion mixing ratios derived from Cassini UVIS data: diamonds 
from the dawn ansa and pluses from the dusk ansa. The solid lines are output from the dual 
hot electron variation model. 



